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ABSTRACT 

Context. The classical flare picture features a beam of electrons, which were accelerated in a site in the corona, hitting the chromo- 
sphere. The electrons are stopped in the dense chromospheric plasma, emitting bremsstrahlung in hard X-rays. The ambient material 
is heated by the deposited energy and expands into the magnetic flare loops, a process termed chromospheric evaporation. In this view 
hard X-ray emission from the chromosphere is succeeded by soft-X-ray emission from the hot plasma in the flare loop, the soft X-ray 
emission being a direct consequence of the impact of the non-thermal particle beam. However, observations of events exist in which a 
pronounced increase in soft X-ray emission is observed minutes before the onset of the hard X-ray emission. Such pre-flare emission 
clearly contradicts the classical flare picture. 

Aims. For the first time, the pre-flare phase of such solar flares is studied in detail. The aim is to understand the early rise phase of 
these events. We want to explain the time evolution of the observed emission by means of alternative energy transport mechanisms 
such as heat conduction. 

Methods. RHESSI events displaying pronounced pre-flare emission were analyzed in imaging and spectroscopy. The time evolution 
of images and full sun spectra was investigated and compared to the theoretical expectations from conduction driven chromospheric 
evaporation. 

Results. The pre-flare phase is characterized by purely thermal emission from a coronal source with increasing emission measure 
and density. After this earliest phase, a small non-thermal tail to higher energies appears in the spectra, becoming more and more 
pronounced. However, images still only display one X-ray source, implying that this non-thermal emission is coronal. The increase of 
emission measure and density indicates that material is added to the coronal region. The most plausible origin is evaporated material 
from the chromosphere. Energy provided by a heat flux is capable of driving chromospheric evaporation. We show that the often 
used classical Spitzer treatment of the conductive flux is not applicable. The conductive flux is saturated. During the preflare-phase, 
the temperature of the coronal source remains constant or increases. Continuous heating in the corona is necessary to explain this 
observation. 

Conclusions. The observations of the pre-flare phase of four solar flares are consistent with chromospheric evaporation driven by a 
saturated heat flux. Additionally, continuous heating in the corona is necessary to sustain the observed temperature. 

Key words. Sun: flares - Sun: X-rays, y-rays - Acceleration of particles 
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The main observational arguments for chromospheric evapora- 
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a\ The question of energy conversion during the impulsive phase ^.^^ ^^^^ ^^.^ ^^^^ ^^^.^^^^ j^^^^ 

O of a solar flare has converged over the course of the past three ^^^^^^ ^^.^^^ 3XRs during flares is one or two or- 

> , decades on a picture featuring an acceleration site in the corona, ^^^^ magnitude greater than is generally measured in die 

where particles accelerated to high energies then precipitate .^^ ^^^^^^ ^ Aschwanden & Acton 2001; Krucker^&U^ 

along magnetic field lines to the chromosphere. In the chromo- ^ ^^.^^^ ^ ^^^^^^ ^j^^ additional material and (2) the 

^ sphere they lose their energy in the dense plasma, leading to ^^^^^^^ upflowing plas ma detected via blueshifts in high 

chromospheric heating, as weU as the characteristic hard X-ray ^ .^^^.^ ti-^in lines (lAntonucci et al.lll982h . in particu- 

footpoint emission. In addition to radiating and conducting its ^^^^^ observations made with imaging spectrometers such 

excess energy away, flie heated chromospheric plasma finds a ^j^^ ^^^.^^^^ Diagnostic Spectrometer (CDS, Hai-risorliEiD 

new equihbrium by expanding up the loop, in a process termed igg gj) the Solar and Heliospheric Observatory (SOHO, 

chromospheric evaporation . The observational evidence for a pomin^o et aLliggS *). Such flows have been observed' in the flare 

relationship between flare heating and non-thermal electi'ons in 41,;^^ j^ase by Miflisan et al. (2006a b) and in the grad- 

flie impulsive phase is reasonably strong Jn this pap er we use ^^^^ by &kowska et^ ([1999^. Evaporative flows have 

m^ging and spectroscopy from the RHESSI satellite ( |LmetaI. | ^^^^ described as being either 'gentle' or 'explosive' , with the 

I2Q02I), as well as information from the GOES satellites to inves- j^^^^^. ^ccuiTing when the rate of chromospheric plasma heat- 

tigate flare heating before the flai-e impulsive phase, appai-endy ■ ^^^^j ^^^^^^^ ^^^^ ^^.^^ ^^^j ^ .^^iation, 

in the absence of non-thermal electrons conduction or expansion, primarily determined by the ratio be- 

Chromospheric evaporation is usuaUy proposed as the cause ^^^^^ heating time scale and the hydrodynamic expansion 

of flare extreme yltr-wiolet i:^^^^^ ^"^ soft X-ray (SXR) emis- ^^^^^^^^^ (iFisher et al.lll985,) . The relatively slow heating rate 
sion (though see IFeIdmanlll990t lActon et alJI 19921 who dispute 



by conduction from a heated corona makes it a lik ely source 
of gentle evaporation. lAntiochos & Sturrockl (Il978h developed 
an analytic model of conductive driven evaporation in the de- 
cay phase of flares, though report that the 'evaporation' repre- 
sented the hydrodynamic redistribution of material already at 
coronal temperatures, rather than chro mospheric materia l being 
heated and expanding into the corona. lKarpen&Devorel(ll987l) 
carried out numerical simulations exploring the effects of non- 
local and saturated heat flux in addition to classical Spitzer heat 
flux. Observational evidence for such "gentle" evaporation was 
found by Za rro & Lemen (JL988) in observations from the Solar 
Maximum Mission, and the flows observe d by iMilligan et al.l 
(12006a') and' Czaykowska et al.l (1 1 99911200 ih were also found to 
be consistent with conductive evaporation, due to both their rel- 
atively low speeds and absence of significant HXR radiation. 

In the electron beam-driven evaporation model, the thermal 
X-ray emission is a direct effect of the energy deposition in the 
chromosphere by the electron beam, so that the time-integrated 
non-thermal hard-X-ray (HXR) flux in a given energy range is 
prop ortional to the SXR flux. This was first proposed by lNeuperti 
( Il968i) and since then, the so-called 'N eupert effect' has been 
studied at length by var i ous authors (e . g.lDennis & Zarrolll993t 
iMcTiernan etaP Il999t IVeronig et all |2005'). Overall the ex- 
pected proportionality between the time-derivative of the SXRs 
and the instantaneous HXRs is supported observationally, par- 
ticularly for impulsive flares, though there is a significant scatter 
interpreted as due to energy losses from low-temperature radia- 
tion, conduction and mass motion. There are also some signif- 
icant e xceptions from this rule - for example, iMcTiernan et al.l 
(1 19991) find that the Neupert relationship manifested in about 
ha lf of the 33 flares the y studied. A similar result was found 
by lVeronig et al.l (l2002ah . This suggests additional energy input, 
not related to non-thermal electrons, in the non-Neupert flares. 
Various authors have also observed impulsive SXR footpoint 
emission tracking th e HXR intensity, rather than its integral (e.g. 
iHudson et al.lll994l) . Finally, pre-flare SXR sources, occurring 
in advance of any HXR emission present a further instance in 
which thermal emission may also be completely unrelated to 
non-thermal electrons. This paper concerns the observation and 
interpretation of such sources. 

It has been known for some time that flare-related activity 
commences prior to the flare impulsive phase. There is a dis- 
tinction to be drawn between pre-flare activity, which refers to 
the very earliest stages of the flare before the impulsive phase 
radiation is detectable, and 'flare precursor' events, which are 
small-scale brightenings in UV to HXR wavelengths happen- 
ing some tens of minutes before the flare. The first to use the 
term "pre-flare" and c omplete a statistical study on flare pre- 
cursors in X-rays were 'Bumba & Kfivskyj (Il959l) . who inferred 
the X-ray behavior via iono spheric disturb ances. With the ar- 
rival of the Yohkoh satellite dOgawara et al.lll99l1) . the nature of 
the pre-flare phase was studied in more detail, particularly the 
relationship between the non-thermal and thermal emissions. It 
was noted relatively early that high temperature thermal sources 
were present in th e corona substan tially before the impulsive 
phase flare onset ( Acton et al.ll 19921) . a result which is clearly at 
odds with the assumption that electron beams drive evaporation. 
Comprehensive statistics including an analysis of the spatial re- 
lations between pre-flares and the consequen t flares h a s been 
conducted by Farnfk et al.. (L996) and Farn fk & Savyl ( |1998| ) 
based on observations wi th the Yohkoh Soft X-ray Telescope 
(SXT. iTsuneta et a"nil99ll) . Their definition of pre-flare activity 
is a rise in the Yohkoh emission above background, five minutes 
to an hour before the main peak. In the spatially-resolved SXT 



images it is clear that in several of their sample there is substan- 
tial coronal soft X-ray emission o ccurring several min utes before 
the start of the impuls ive phase ( | Farnjk & Savvll998l) . In a study 
of 10 Yohkoh flares, lAlexander et al.l (1 19981) also identify line 
broadening observed with Yohkoh Bragg Crystal Spectrometer 
(Culhane et al. 199JJ several mi n utes b efore the onset of the 
impulsive phase, and iHarra et aP (1200 ll) show that in one such 
event the temperature, intensity, and non-thermal line width are 
changing significantly, well before the start of the hard X-rays. 
Taken together, the Yohkoh observations provide strong evidence 
for pre-flare coronal activity unrelated to evaporation driven by 
electron beams. In recent years, RHESSI observations further 
confirmed the existence and relative r elative frequency of pre- 
flare activity as shown eg. in lVeronig^ t al. (2002b). 

The Ramaty H igh Energy Solar Spectroscopic Imager 
(RHESSI. lLin et al.ii2002) offers new opportunities to study pre- 
flare sources in detail. RHESSI is a HXR telescope, with spectral 
coverage extending down in energy to the thermal range at ~ 3- 
20 keV. Its high spectral resolution (of about 1 keV at energies up 
to 100 keV) permits spectroscopic diagnostics not available pre- 
viously, coupled with high spatial resolution (as low as 2.3", de- 
pending on the flare) in user-defined energy channels. This offers 
new temperature and emission measure diagnostics for very hot 
plasmas (around 10 MK and above) as well as for non-thermal 
electrons, if present. The RHESSI observations presented here 
display an increase of low energy X-ray emission, which can be 
fitted with a thermal spectrum, up to minutes before the rise of 
the HXR emission. 

The morphology of these events shows only one source, 
visible at energies in the thermal regime, while the footpoints 
only appear at the onset of the HXR emission. We will argue 
that this indicates that another mechanism of energy transport to 
the chromosphere may be important, such as heat conduction. 

In this study, we present the first comprehensive RHESSI-led 
study of the pre-flare phase of solar flares. The paper is struc- 
tured in the following way: Section |2] describes the flare selec- 
tion process. In Sect.|3]the time evolution of the selected flares is 
analyzed in spectra and images. A theoretical model to explain 
the observed time evolution is presented in Sect. HI followed by 
a Discussion and Conclusions. 



2. Flare selection 

In this study, events with pre-flare activity will be analyzed. 
The main pre-requisite of such a flare is increasing emission in 
SXR well before the onset of the HXR emission. Therefore, only 
events the beginning of which was fully observed were consid- 
ered. The selection criteria are summarized here: 

- Offset of more than 700" from disk center to limit projection 
effects in the study of the flare morphology; 

- GOES class larger than Ml to ensure good count rates for 
imaging; 

- The beginning of the event had to be observed with RHESSI; 

- SXR & GOES emission had to start increasing at least 1 min 
before the HXR emission; 

- Events with a simple morphology involving one SXR source 
at the onset and two HXR sources appearing later in the 
event. This indicates the time when the acceleration process 
becomes dominant; 

- High enough detector livetime (> 90 %) for reliable spec- 
troscopy and imaging. 
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Fig. 1. RHESSI corrected count rate lightcurves in 6-12 keV and 25-50 keV during the pre-flare and early impulsive phase of the 
events. Missing data in the RHESSI lightcurves indicate the time interval that was omitted in the analysis due to attenuator state 
change. The dotted vertical lines indicate the time interval of the primary analysis, the dashed line gives the start time of the image 
in which footpoints are observed first. GOES lightcurve (dash-dotted). 



Table 1. Eventlist. The attenuator state is for all events at the very beginning, then changes to 1 as the flux increases. The attenuator 
times give the rough time interval that must be avoided because of attenuator motion. 



Date Analyzed time GOES class xy-position [arcsec] coronal source volume [cm^] Attenuator Attenuator times 

03-Aug-2002 18:59-19:02 XO 912/-271 4.9-10^" 0/1 19:02:20/19:03:10 

13-NOV-2003 04:55:40-04:57:12 M1.7 -977/23 0/1 04:57:20/04:58:20 

06-Jan-2004 06:13-06:16 M5.9 -992/88 2.8-102' 0/1 06:18/06:19 

21-May-2004 23:37:30-23:39 M2.6 -757/- 157 1 0/1 23:40:30/23:41:30 



Four events were selected as best suited for this kind of study. 
Their key parameters are listed in Table[T] The RHESSI attenua- 
tors were out during the earliest phase of all events with the thin 
attenuator moving in as the count rates increased, typically after 
2-3 minutes. 

3. Time evolution of events 

The time evolution of the SXR emission in GOES and RHESSI 
6-12 keV along with the HXR emission (RHESSI 25-50 keV) 
is shown in Fig.[T] The missing data in the RHESSI lightcurves 
indicate the time intervals omitted due to attenuator state 
change. The attenuator state was for all events at the onset of 
the flare. 



pulsive phase were made, using grids 3 to 8. The time interval 
of each image was 30s in order to get high enough count rates 
throughout the pre-flare phase. Full-sun photon spectra were 
fitted in 30 s time intervals, corresponding to the image time 
intervals. The fitting model consisted of a thermal component 
fitted from 6 keV and, if possible, a non-thermal component. 
Spectra from the onset of the flares as well as the start of 
the non-thermal emission are shown in Fig. |4] The RHESSI 
spectrum at the lowest energies (3-6 keV) could be fitted with 
an additional thermal component with temperature and emission 
measure similar to the values that can be derived from GOES 
measurements (not shown in the figure). 



3.1. Data analysis 

We studied the time evolution of the pre-flare phase both in 
images and spectra. CLEAN and Pixon images dHurford et al.l 
l2002t iMetcalf et aP [19961) during flie pre-flare and early im- 



3.2. Spatial evolution 

Figure 12] shows CLEAN images of the events taken at the begin- 
ning, at the time when a non-thermal component was first fitted 
and at the time when the footpoints first appear For each event, 
3 energy-bands (6-12, 12-25, 25-50 keV) are shown. 
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Fig. 2. Contour plots from CLEAN images (60%, 75% and 90% of maximum flare emission) of all events taken at the first fit interval 
of the rise phase (left), the time interval when a non-thermal component could first be fitted {middle) and at the first appearance of 
the footpoints (right). Those times correspond to the times of the spectra shown in Fig. HI The energy bands are 6-12 keV (top), 
12-25 keV (middle) and 25-50 keV (bottom). 



All events start with a single source visible only at the lowest 
energies in the range 6-12 keV (upper left image of the 9-image 
panel per flare in Fig. [21 This is interpreted as a source at the top 
of a loop (coronal source). After some minutes two additional 
sources appear at higher energies (25- 50 keV) which are inter- 
preted as chromospheric footpoints (lower right image of each 
flare in Fig.|2]i. In three events the position of the first appearing 
source is clearly displaced from the footpoint position, imply- 
ing that there is an actual loop geometry with a coronal source 
on top and footpoints in the chromosphere. In the fourth event, 
the separation is not as clear but from the spectral evolution (see 
next Section) one can still assume a coronal-source-footpoints- 
geometry. 



3.2.1 . Movement of coronal sources 

We analyzed the time evolution of the source position by mea- 
suring the centroid position of the 50% contour in CLEAN and 
Pixon images in the energy range 6-12 keV. The position from 
both imaging methods agree within the uncertainties. Figure |3] 
displays the 70% contour from CLEAN images with positions 
as found from CLEAN. Colors go from black to light (time of 
first appearance of footpoints). The positions shown indicate the 
position of the coronal source from the beginning until the time 
of the first appearance of footpoints, in the same time steps as 
used in spectroscopy (30s). 

In all events, the sources move to some extent during the 
pre-flare phase. The source in the event of 03-Aug-2002 is sta- 
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Fig. 3. Position of the coronal source in time in the energy range 
6-12 keV. Time evolution is indicated by 70% intensity contours 
in colors ranging from black to light with the first appearance of 
the footpoints (25-50 keV) given as dotted areas in the appro- 
priate color The centroid position is marked with a cross in the 
corresponding color 



ble within the position uncertainties during the pre-flare phase, 
but appears at a displaced location in the last image (at the first 
appearance of the footpoints). Due to the attenuator state change 
just before the last image, an interval of about 1 minute is miss- 
ing. Therefore, one cannot say whether the source moves contin- 
uously or not. The same observation is made in the 13-Nov-2003 
event. Here, the sources observed during the pre-flare phase and 
at the beginning of the impulsive phase c ould possibly i ndicate 
two separate loop systems as suggested bv lLiu et"an (2006). The 
source position of the 6-Jan-2004 event remains constant over 
time. The 21-May-2004 event is somewhat peculiar, display- 
ing pronounced, continuous source motion. The displacement 
of the position is clearly larger than the uncertainties. From the 
direction of the movement, this could be another case of alti- 
tude decrease as has been found in othe r RHESSI events (eg. 
LSui & Holmanll2003HVeronig et alJl2006h 

3.3. Spectral evolution 

All four events display a purely thermal spectrum in the first of 
the analyzed time intervals. This spectroscopic finding is fur- 
ther supported by GOES SXI images (not shown here). These 
are available for three out of four events and show a soft X- 
ray source minutes to hours before the RHESSI observations. In 
RHESSI observations the apparently purely thermal phase lasts 
one to two minutes, after which a small tail to higher energies 
becomes visible (Fig.|4]l. At this stage, there are as yet no foot- 
points are observed, indicating that the non-thermal emission is 
coronal. The tail appears independent of albedo correction and 
as the events are all near the limb, the influence of albedo on 
the total flux is minimal. Further, pile-up is not large enough at 
this stage to account for the emission. It is therefore safe to as- 
sume that the tail is real. It can be fitted with a power-law or 
with a second, very hot thermal component. During the impul- 
sive phase (after attenuator state change) the HXR component 



becomes very pronounced and clearly distinguishable. Figure |4] 
shows spectra of the selected events, one taken at flare onset and 
one when the non-thermal tail first appears. 

3.4. Time evolution of flare parameters 

The time evolution of the fit parameters, namely temperature, 
emission measure and spectral index is shown in Fig.|5] The fig- 
ure displays panels for each event. The top panel of each event 
shows RHESSI lightcurves in the 6-12 and 25-50 keV energy 
band. The cross indicates the time of first appearance, and the 
hardness, of the non-thermal component. The second panel dis- 
plays the emission measures as determined from RHESSI spec- 
tra and GOES. In the third panel, temperatures are presented the 
same way as the emission measures. The fourth panel shows 
the electron density estimates. The critical parameter in com- 
puting densities is the flare volume. RHESSI images allow an 
estimate of the volume by measuring the source area A and ap- 
proximating the volume as V = A'^^^. Measuring the area from 
RHESSI images is quite imprecise with large uncertainties. We 
determined the 50 % contours in CLEAN images (de-convolving 
the CLEAN beam). The area is constant over the pre-flare evo- 
lution and we made an estimate of the quantitative value of the 
volume by taking the time average of the measured values. The 
average electron density in the coronal source is then given as 

^ ^JEM/V, (1) 
assuming a volume filling factor of 1 . 

4. A theoretical model to explain the observations 

In Fig.|5]we present the time evolution of flare parameters such 
as temperature, emission measure and density, as well as im- 
ages. The observations reveal increasing emission measure and 
densities over the course of the pre-flare as well as generally in- 
creasing temperatures. Assuming a constant, or even increasing, 
source volume, the density increase can only be attributed to ad- 
ditional material which is added to the coronal source region. 
The most likely scenario causing such a density increase is chro- 
mospheric evaporation. In this scenario, chromospheric plasma 
(including the transition region) is heated and expands upward 
the magnetic loops. There are two possible main heating mech- 
anisms for the chromosphere; non-thermal electron beams and 
thermal conduction. In the former case, a beam of high energetic 
electrons impinges on the chromosphere, where the particles are 
completely stopped in the dense target. As a consequence of this 
energy deposition, the temperature rises. In the latter case, en- 
ergy is transported from the hot coronal source to the cooler 
chromosphere by thermal conduction. We are going to discuss 
the thermal conduction scenario further in the following subsec- 
tion. Reasons why thermal conduction is favored over particle 
beams are presented in the discussion. 

4. 1 . Theory of thermal conduction 

Our model will explain chromospheric evaporation by heating 
due to energy input from heat conduction. The coronal source is 
much hotter than the chromospheric plasma, therefore a temper- 
ature gradient develops, driving a heat flux downward along the 
magnetic loop. 

The conductive heat flux is given as 

Fcond = KOT'I^^. (2) 
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Fig. 4. Spectra with model fits for 2 time intervals. Top: First analyzed time interval. Bottom: First interval during which a non- 
thermal component could be fitted. 



In classical thermal conduction theory foUowing'Sp itzej (Il965h . 
the conductive coefficient is kq = 10"^ erg cm^'s^'K"^**^. For 
a flare of half loop length Lioop and coronal source temperature 
Tcs, this may be approximated by 



cond 



10 



7,7/2 

-6 ^ cs 
-'loop 



[erg cm ^s ']. 



(3) 



However, the maximum heat flux a plasma can carry is lim- 
ited by a fraction of the thermal energy-flux F,f,erm - nekTeV,h, 
where y,/, is the thermal electron velocity. For sufficiently large 
temperature gradients, such as occur in solar fla res, the heat flux 
is exp ected to reach this limit and saturate. Grav & Kilkennv 
(119801) showed that there are two regimes of flux saturation. 
A non-classical treatment with local flux limiting is nec essary 
if the electron mean free path Aemj - 5.21 ■ IQ^T^/tie (iBenj 
I2OO2I) is larger than only 0.12% of the temperature scale length 
Lii, . If the electron mean free path even exceeds the temperature 
scale length, the heat flux becomes non-local in the sense that it 
depends on the global density and temperature structure of the 
plasma. Figure |6] illustrates the boundary conditions for locally 
limited and non-local heat flux in the density-temperature plane 
of typical solar flare values. The locally limited regime is de- 
scribed bv lCampbelll (11984 as a continuous transition from the 
purely classical treatment to ffie non-local regime. Depending 
on the ratio 'R - Aemf/L,/,, a reduction factor £) < 1 is applied 
to the classical heat flux, resu l ting in an eff'ective heat flux of 
Fred - Q{fi) ■ Fcond- ICampbell (1 19841) published the values of g 
for different H and atomic number Z. Using the values for Z=l, 
we fitted a function of the form 

p = A ■ e-*(*+'^'' (4) 

to the values published bv lCampbelll (|1984|) . The parameters are 
X = InC??), A=1.01, b=0.05, c=6.63. Using those values, the lo- 



cally limited, reduced heat flux can be computed for any H from 
Eq. dUi. 

In the non-local regime where the electrons are freely 
streaming and complete saturation is reached, the flux can be 
expressed as; 



Fja, = 0.53n(,mfy^,j erg cm ^s ' (5) 
(iKarpen & Devordl 19871) . 



4.2. Observation of thermal conduction 

In studies of ffiermal conduction in flare loops, a classical sce- 
nario involving Spitzer conduction is often assumed. Are the 
conditions in the observed pre-flares consistent with classical 
conduction or is the flux limited? The densities and tempera- 
tures as derived from GOES and RHESSI observations are dis- 
played in Fig.|6| Values measured with GOES are given as dots, 
RHESSI measurements are marked as triangles. GOES typically 
yields larger emission measures and smaller temperatures than 
RHESSI. Assuming the same emission volume, this results in 
higher densities for GOES. The figure reveals that the condi- 
tions during the observed times of the pre-flare cannot sustain a 
classical heat flux, the heat flux is limited. The resulting effec- 
tive flux was computed from the observed values. Temperatures 
and densities yield the electron mean free path A,„fp. The tem- 
perature scale length L,/, is given as L,i, = T/AT x Lioop, where 
Ar was approximated as AT ^ Tcs/Lioop. Using an estimated 
loop half-length of 10^ cm we computed the ratio fi. The cor- 
rection factor and therefore the effective heat flux was computed 
from Eq. (|4]i. The results are displayed in Fig.|2l again for both, 
GOES and RHESSI values. 
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Fig. 5. Time evolution of the fit parameters for all 4 events. Top: RHESSI lightcurves in the 6-12 and 25-50 keV energy band. The 
cross indicates first appearance and hardness of the non-thermal component. Second: Emission measure as computed from RHESSI 
{crosses and GOES dash-dotted). Third: Temperatures, same as for emission measures. Fourth: Electron density in the coronal 
source, same as for emission measures. 



4.3. Chromospheric evaporation 

We assume that the heat flux computed in the section above 
is deposited in the transition region and chromosphere, heating 
the chromospheric material. The hot material will then expand 
upward the magnetic loops. This chromospheric evaporation is 
usually divided into two types; gentle evaporation and explosive 
evaporation. The determining factors distinguishing the two are 
the ratios between heating, radiative cooling and expansion rate 
in the chromosphere. If the heating rate is much smaller than the 
radiative loss rate, all the deposited energy is radiated away. If 
the heating rate is only marginally larger than the radiative loss 
rate, the temperature and pressure rise slowly and the material 
evaporates gently. If the deposited energy rate exceeds a certain 
threshold value, the chromospheric plasma is heated to coronal 
temperatures quickly, a large overpressure develops and explo- 



sive expansion of plasma is the consequence. iFisher et alJ(ll985h 
derived a value for the critical heating flux of 10'° erg cm^^s"' 
for typical coronal and chromospheric conditions. It is gener- 
ally assumed that explosive evaporation is associated with the 
large energy input provided by beams of fast electrons hitting the 
chromosphere, while gentle evaporation is normally attributed to 
heating by thermal conduction. 

4.3.1. Conduction driven cliromosplieric evaporation 

Can the observed emission measure and density increase in the 
presented events be explained with chromospheric evaporation 
driven by the heat flux presented in Sect I4.2P The total energy 
flux in the evaporated plasma is given as 

F:;r = 3nr'kr^-'\„ = 3myrvlv,., [ergcm-2s-i], (6) 




^q8 ^qIO ^Qi2 

Electron density [cm^^] 



Fig. 6. Temperature vs. density space of typical flare loop values, 
indicating the boundaries and regimes for classical thermal con- 
duction, flux limited and non-local conduction. The density and 
temperature values of the observed events are given as dots for 
GOES and triangles for RHESSI observations. Different shades 
stand for the different events. 



We use the evolution of the density in the coronal source to 



estimate v„p. 
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The density of the upflowing plasma cannot be determined from 
observations. Assuming n^"'' - n^, the density in the coronal 
source, the upflow velocities derived from Eq. (O range from 
20 km s"' to 150 km s"'. This results in an evaporative flux of 
the order of 10^ erg cm^^s ' . Further, the evaporative expansion 
requires an energy of 

Ee,p = J pad\ (9) 

where p is the total pressure. Assuming isoth ermal expansion 
the ex pansion energy flux can be computed (see lBenz & Kruckeil 
l2002i) : 

F exp — ^f^ef^e^th^up MLiooplU) (10) 

The expansion energy flux is larger than the thermal energy flux 
(Eq.|6]l by a factor of 2/3ln(Lioop/li) which amounts to about 5. 

The total energy flux for evaporation is therefore in the range 
of 5 • 10^ erg cm^^s ', which is well below the values obtained 
from the data (Fig.|7]i. The rest of the conducted energy heats the 
transition region and upper chromosphere to sub-coronal tem- 
peratures and is radiated away. 




Fig. 7. Effective conductive flux derived from RHESSI (crosses 
on solid line) and GOES (solid line). The dottet line marks 
the threshold ener gy flux to trigger explosive evaporation (from 
iFisher et alJflQSSh . The thin dashed lines mark the level of un- 
certainty (see|5l). 



where y„p is the velocity of the upward moving plasma, v,;, the 
thermal velocity and the density of the upflowing plasma. 



5. Discussion 

5. 1 . Electron beams in the pre-flare phase 

In the above section, we demonstrated that conduction driven 
chromospheric evaporation can account for observations of pre- 
flares. Can electron beams below the detection limit do the 
same? Let us estimate the maximum expected non-thermal elec- 
tron flux and the expected chromospheric heating rate. We esti- 
mate the electron beam strength starting from the observed coro- 
nal source spectrum. The first time-interval when non-thermal 
emission from the corona was observed was used as an esti- 
mate of the upper limit to the electron flux, assuming a thin 
target in the coronal source. If a beam constituting of this elec- 
tron flux hits the chromosphere, a thick target footpoint spec- 
trum Ffp - Aii is exp ected. From observation s of footpoints 
and coronal sources (eg. iBattaglia & Benzll2006l) . we know that 
a flux Ffp as found from this estimate would be well observed 
in images and spectra which is not the case. Taking the obser- 
vations of Battaglia & Be nz (2006) as a reference, we estimate 
that the effective electron flux that reaches the footpoints has to 
be about a factor 100 smaller than expected from the coronal 
source observations. For the computation of the beam heating 
rate, we therefore use an assumed footpoint photon flux spec- 
trum F'/J = 0.01 ■ Ffp. 

The heati ng r ate of a beam of ele ctrons was derived by 
iBrownl (Il973h and lLin & Hudsoi] (Il976h . The rate of energy in- 
put per unit volume at column density and time t is given as 

lB(N,t) = 10\(N, t)[x(N,t) + 0.55] 

X -C(r)(l.l ■ 10"'^A^)"'''+'>^^ [ergcm^^] (11) 
a 



where n{N, t) is the hydrogen density, x{N, t) the fractional ion- 
ization, a the cross-section area, A the normalization of the ob- 
served photon spectrum, y the spectral index of the observed 
photon spectrum and C{y) a function of the y6-function and y. 
The parameters of the expected footpoint spectrum were used as 
input for A and y in Eq. ( fTTT l. As an estimate of the column den- 
sity, the stopping d epth of a 35 keV electron was used a ccording 
to A' = lO^^ii? dTandberg-Hanssen & Emslidll988l) . The re- 
sulti ng density n(A^ , t) and fractional ionization x{N, f) was taken 
from lBrownl d 1 973. Tablel). The cross-section area was taken as 
lO'^ cm^, which corresponds to typical footpoint areas as mea- 
sure d by RHESSI. Finally, the factor C{y) was taken from Fig. 
2 in iLin & HudsonI (Il976h . Inserting the above values gives a 
total energy input of the order of 3-4 erg cm^^s From the 
expression of the column density = nl, where / is the path 
length, 1 can be determined and by multiplying Eq. (fTTT) with 
/ we get an energy input of 5-7-10^ erg cm^^s"', which is two 
orders of magnitude smaller than the energy input by thermal 
conduction. Even if the effective footpoint area was an order of 
magnitude sm aller, as often found in UV or Ha observations (eg. 
iTemmer et alji200 7). the heating rate by the beam would still be 
an order of magnitude smaller than the energy input by thermal 
conduction. Therefore, if there is any chromospheric heating by 
electron beams in the pre-flare phase, its contribution is minimal 
compared to the energy input by thermal conduction. 

5.2. Influence of assumptions 

In all computations throughout this work, the values measured 
in the coronal source were used also for the loop legs. For some 
parameters, e.g. the electron mean free path /!/„//„ those values 
may not be accurate. For a complete treatment, the temperature 
and density in the loop is needed. However, those values can- 
not be determined observationally. If the loop was at the same 
temperature and density as the coronal source, it would be vis- 
ible in X-ray images. As no loop emission is observed in the 
images, the actual loop has to be either less dense or cooler than 
the coronal source. Thus the coronal source values represent an 
upper limit to the true conditions in the loop. Smaller temper- 
ature could result in an order of magnitude smaller fi. On the 
other hand, a smaller density could cause an order of magnitude 
larger ratio fi of electron mean free path to temperature scale 
length. The correction factor g (Eq. HI would change by 0.5 to 
2, respectively and so would the conductive flux. The resulting 
upper and lower limits of the conductive flux are indicated by 
dashed lines in Fig. [7] Even if the conductive flux differed by a 
factor of 2, this would not change the basic interpretation. 

The evaporated flux F"p''' ~ iieTvup ~ T only depends on 
the temperature. It is therefore expected to be at most about a 
factor of 4 smaller. However, the upflow velocities derived from 
Eq. ^ are proportional to 1 /n^, and could therefore be up to an 
order of magnitude larger for smaller densities. Eq. (O may then 
give velocities higher than the local sound speed which would 
be unphysical. 

5.3. Mode of conduction 

Our results indicate that the conductive heat flux during the pre- 
flare phase is non-classical. The temperature gradient between 
the coronal source region and the chromosphere is so steep that 
the classical Spitzer conduction is not valid. The conductive flux 
is locally limited for temperature and density values determined 
from GOES and may be fully saturated (i.e. the electrons are 



freely streaming) in the case of the RHESSI measurements. Both 
situations lead to a reduced heat flux compared to the classical 
Spitzer conductivity. This finding does not change directly the 
conclusions drawn from the work presented here. However, it re- 
quires that care be taken when computing total flare energy bud- 
gets as the energy loss of the coronal source due to conduction 
may be smaller than anticipated from the assumption of classical 
Spitzer conductivity. 

5.4. Acceleration vs. heating 

Despite the energy loss in the corona due to the heat flux, in- 
creasing temperatures during the pre-flare phase are observed. 
Continuous heating in the coronal source is therefore necessary 
to sustain and even increase the observed temperature. Possible 
heating mechanisms include heating at the stand-off slow-mode 
shocks in Petschek reconnection models which has been pro- 
posed as an explanati on of the fo rmation of hot thermal sources 
in Yohkoh ff ares (eg. iTsuneta e t a l. 1997), or betatron hea ting 
in a collapsing magnetic trap (eg. iKarlicky & Kosugill2004l) . A 
strong candidate heating mechanism is wave-particle interac- 
tion as described in the transit-time-damping model ( Mill er et all 
1996; Grigis & Benz 2006). Such a model could also explain the 
spectral time evolution of the observations: Depending on the 
escape time and wave density in the acceleration region, parti- 
cles are not accelerated to non-thermal energies but remain at 
a quasi-Maxwellian distribution observationally not distinguish- 
able from a purely thermal plasma at enhanced temperature. For 
increasing wave density or escape time, particles are acceler- 
ated to higher energies, constituting a non-thermal tail which is 
clearly noticeable in the observations. 

However, the three dimensional geometry of the flare site 
has to be considered and a smooth transition from heating to 
acceleration at the same site may not reflect reality in two of 
the presented events. As shown in Fig. [3] the position of the 
coronal source changes over time in the events of 03-Aug-2002 
and 13-NOV-2003. Farniketal. (1996) and Fai-nik & Sav^ 
(il998l) made a statistical study on the position of Yohkoh 
pre-flares and the successive impulsive flares, finding that only 
in 25% of the observed events were the pre-flare emission and 
impulsive flare emission spatially coincident. The other events 
were classified into what the authors termed "distant" and 
"adjacent/overlapping" events. In the latter case, parts of the 
pre-flare 50% maximum intensity emission overlap with parts 
of the main flare emission while in the "distant" events, the 
two emission patterns are spatially separated. If our two flares 
indeed fall into the category of distant events, it will not change 
the interpretation of the pre-flare phase, but care should be taken 
when interpreting the time-evolution of the spectrum beyond 
the purely thermal phase. 



6. Conclusions 

The observed increase in thermal emission and SXR emission 
measure up to minutes before the start of the impulsive flare 
phase in the four events presented here can be explained by 
chromospheric evaporation driven by thermal conduction. To 
compensate for the heat loss, continuous heating in the coronal 
source is needed. In this early phase, acceleration of particles to 
non-thermal energies is minimal or non-existent. A scenario fol- 
lowing stochastic acceleration is conceivable in which the coro- 
nal source region is only heated in the early flare stages. As soon 
as the acceleration mechanism becomes efficient enough, a tail 



of non-thermal particles is produced, first visible in the coronal 
source and eventually in the footpoints. 
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